Computation of Gauss-kronrod Quadrature Rules with Non-positive Weights
نویسنده
چکیده
Recently Laurie presented a fast algorithm for the computation of (2n + 1)-point Gauss-Kronrod quadrature rules with real nodes and positive weights. We describe modifications of this algorithm that allow the computation of Gauss-Kronrod quadrature rules with complex conjugate nodes and weights or with real nodes and positive and negative weights.
منابع مشابه
Computation of Gauss-Kronrod quadrature rules
Recently Laurie presented a new algorithm for the computation of (2n+1)-point Gauss-Kronrod quadrature rules with real nodes and positive weights. This algorithm first determines a symmetric tridiagonal matrix of order 2n+ 1 from certain mixed moments, and then computes a partial spectral factorization. We describe a new algorithm that does not require the entries of the tridiagonal matrix to b...
متن کاملOrthogonal Polynomials and Quadrature
Various concepts of orthogonality on the real line are reviewed that arise in connection with quadrature rules. Orthogonality relative to a positive measure and Gauss-type quadrature rules are classical. More recent types of orthogonality include orthogonality relative to a sign-variable measure, which arises in connection with Gauss-Kronrod quadrature, and power (or implicit) orthogonality enc...
متن کاملCalculation of Gauss-Kronrod quadrature rules
The Jacobi matrix of the (2n+1)-point Gauss-Kronrod quadrature rule for a given measure is calculated efficiently by a five-term recurrence relation. The algorithm uses only rational operations and is therefore also useful for obtaining the Jacobi-Kronrod matrix analytically. The nodes and weights can then be computed directly by standard software for Gaussian quadrature formulas.
متن کاملQuadrature rules for rational functions
It is shown how recent ideas on rational Gauss-type quadrature rules can be extended to Gauss-Kronrod, Gauss-Turr an, and Cauchy principal value quadrature rules. Numerical examples illustrate the advantages in accuracy thus achievable. 0. Introduction The idea of constructing quadrature rules that are exact for rational functions with prescribed poles, rather than for polynomials, has received...
متن کاملStable Computation of High Order Gauss Quadrature Rules Using Discretization for Measures in Radiation Transfer
The solution of the radiation transfer equation for the Earth's atmosphere needs to account for the re ectivity of the ground. When using the spherical harmonics method, the solution for this term involves an integral with a particular measure that presents numerical challenges. We are interested in computing a high order Gauss quadrature rule for this measure. We show that the two classical al...
متن کامل